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Second-order power spectra of Cosmic Microwave Background (CMB) anisotropies due to random 
primordial perturbations at the matter dominant stage are studied, based on the relativistic second- 
order theory of perturbations in flat cosmological models and on the second-order formula of CMB 
anisotropies derived by Mollerach and Matarrese. So far the second-order integrated Sachs- Wolfe 
effect has been analyzed using the three-point correlation or bispectrum. In this paper we derive the 
second-order term of power spectra given using the two-point correlation of temperature fluctuations. 

The second-order density perturbations are small, compared with the first-order ones. The second- 
order power spectra of CMB anisotropies, however, are not small at all, compared with the first-order 
power spectra, because at the early stage the first-order integrated Sachs- Wolfe effect is very small 
and the second-order integrated Sachs- Wolfe effect may be dominant over the first-order ones. So 
their characteristic behaviors may be measured through the future precise observation and bring 
useful informations on the structure and evolution of our universe in the future. 

PACS numbers: 98.80.-k, 98.70.Vc, 04.25. Nx 

I. INTRODUCTION 

In most studies of Cosmic Microwave Background (CMB) anisotropies, the comparison between observed and 
theoretical quantities have so far been done, assuming the linear approximation for cosmological perturbations. 
It seems to be successful enough to determine the cosmological parameters [l], 0, H| • The present state of our 
universe is, however, locally complicated and associated with nonlinear behavior on various scales, and so the 
observed quantities of CMB anisotropies may include some small effects caused through nonlinear process by 
various primordial perturbations. So far the nonlinearity in CMB anisotropies has been treated as the second- 
order integrated Sachs- Wolfe effect, which was analyzed using the three-point correlation and bispectrum. 0,0,0 
Recently a general treatment of second-order temperature fluctuations has been systematically studied by 
Bartolo, Matarrese and Riotto@,[l] due to the transfer function which was derived from the Boltzmann equation, 
so as to be applied not only to the integrated Sachs- Wolfe effect at the matter-dominant stage, but also the 
nonlinear effect at the recombination epoch and the primordial nonlinear effect. 

In recent years, on the other hand, we have studied these nonlinear effects of inhomogeneities on CMB 
anisotropies at the matter-dominant stage, based on the relativistic second-order theory of cosmological 
perturbations. 3 and on Mollerach and Matarrese's second-order formula of CMB anisotropies [l(]| . In previous 
papers. [Til fl2l.Tl3l. Hil flU we have studied the second-order effects of a local special large-scale inhomogeneity 
on CMB anisotropies, paying attention to the interaction between it and primordial random perturbations. In 
this relativistic theory we used the first-order and second-order perturbations in pressureless matter which were 
derived in compact and analytic forms. In the present paper we study the nonlinear effect of only primordial 
random perturbations on CMB anisotropies and derive the second-order power spectra, which is given as the 
second term in the two-point correlation of CMB anisotropies. This is a nonlinear correction to the first-order 
power spectra and different from the three-point correlation and bispectrum which have been used in the above 
other works to investigate uniquely the non-Gaussianity in the perturbations generated in various stages. 

The second-order random density perturbations themselves are small, compared with the first-order ones. 
In the A-dominated model, however, the second-order power spectra of CMB anisotropies are not small at all, 
compared with the first-order power spectra, because at the early stage the first-order integrated Sachs- Wolfe 
effect (ISW) is very small and the second-order ISW may be dominant over the first-order one. Especially 
in the local void model with A = o[H, [TH, H3, HH, the exterior region is described by use of the Einstein- 
de Sitter model, in which the first-order ISW vanishes and the second-order ISW is indispensable. So their 
characteristic behaviors may be measured through the future precise observation and bring useful informations 
on the structure and evolution of our universe in the future. 

In Sec. II, we show the second-order perturbations in general flat cosmological models and the corresponding 
CMB anisotropies. In Sec. Ill, we derive the expressions for the second-order power spectra of CMB anisotropies, 
due to primordial random perturbations. Concluding remarks follow in Section IV. 
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II. SECOND-ORDER METRIC PERTURBATIONS AND TEMPERATURE ANISOTROPIES 



First we review the background spacetime and the perturbations which were derived in the previous paper [ll|. 
The background flat model with dust matter is expressed as 



ds 2 — g liV dx^dx v — a 2 {rf)[—drj 2 + Sijdx t dx-'], 



(2.1) 



where the Greek and Latin letters denote 0, 1, 2, 3 and 1, 2, 3, respectively, and (= <5* = 5^) are the Kronecker 
delta. The conformal time rj(= x°) is related to the cosmic time t by dt — a{rf)drj. The matter density p and 
the scale factor a have the relations 



pa 2 = 3(a'/a) 2 — Aa 2 , and pa 3 = po, 



(2.2) 



where a prime denotes d/drj, A is the cosmological constant, po is an integration constant and the units 
8ttG = c = 1 are used. 

The first-order and second-order metric perturbations Sig^v(= h^ v ) and 82 9^v(= respectively, were 

derived explicitly by imposing the synchronous coordinate condition: 



hoo =h oi = and £ Q0 = £ Ql = 0. 
Here we show their expressions only in the growing mode: 

hi = p^Fij, 

i\ = P(r t )L\+P 2 ( J1 )Ml+Q(r t )N\>+Cl, 



(2.3) 



(2.4) 



where F is an arbitrary potential function of spatial coordinates a; 1 , a; 2 and x 3 , h{ = S^hn, N'J = S^N^n 
F^j means d 2 Fjdx % dx^ and P{j]) and Q{rf) satisfy 



a 

Q" + —Q' = -\p-\{p') 2 

a L 2 



(2.5) 



The three-dimensional covariant derivative \i are defined in the space with metric dl 2 = 8ijdx' l dx : > and their 
suffixes are raised and lowered by use of 5%j . The functions L\ and M\ are defined by 



L 



3 _ 



:'./ ',/•:, - 2F ■ F :ij + -SijFjFj 



M j = — { 
28 I 



19F. iZ F. 



and N is defined by 



AN 



12F«AF-3& 



FmFm 



(AF) 2 ]} 



28 



(AF) 2 - F M F t 



hi 



The last term C\ satisfies the wave equation 

3 



-(P/a) 2 Gj + ±[p-l(pi) 2 



(2.6) 



(2.7) 



(2.8) 



where G\ and G\ are second-order traceless and transverse functions of spatial coordinates, and the operator 
□ is defined by 



□ = gThuv = ~a- 2 (d 2 /d V 2 + ^d/d V - A) 0, 



(2.9) 



where ; denotes the four-dimensional covariant derivative. So C\ represents the second-order gravitational waves 
caused by the first-order density perturbations. 
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The velocity perturbations Si and S2 u M vanish, i.e. Si u° = Si u 1 — and S2 u° — S2 u z = 0, and the 
density perturbations are 

SP/P = ^(-P'-lW, 



pa 2 \ a 



Sp/P = 1 r^{l(l--P')(3F >l F i i + 8FAF) + l-P[(AF) 2 +F <kl F <kl ] 
2 2pa z 1 2 a 2 



1 



(P'f - l-Q'] [(AF) 2 - F kl F kl ] - l-PP'[4F kl F kl + 3(AF) 2 ]}. (2.10) 
7 a J 7 a J 



Next let us consider the CMB temperature T = T(°)(l + ST/T), in which is the background temperature 
and ST/T(— SiT/T + S2T/T) is the perturbations. The present temperature To^ is related to the emitted 

background temperature Tj ^ at the recombination epoch by — (1 + z e )To°\ the temperature perturbation 
t = (ST/T) e at the recombination epoch is determined by the physical state before that epoch, and the present 
temperature perturbations (ST/T) D is related to (5T/T) e by the gravitational perturbations along the light ray 
from the recombination epoch to the present epoch. The light ray in the unperturbed state is described using the 
background wave vector (= dx^/dX), where A is the affine parameter, and its component is fc^ ^' = (1, — e 1 ), 
and the ray is given by x^ ^ = [A, (Ao — A)e 4 ], where e l is the directional unit vector. Here and in the following 
the suffixes o and e for A, rj and r denote the present (observed) epoch and the recombination (emitted) epoch, 
respectively. 

The first-order temperature perturbation is 

ST/T = t+ \ ( ' dAP'^F^eV. (2.11) 



1 2 - 

Using the relations dP/dX = P' and dF/dX = —F^e 1 , this equation (|2.11[) is expressed as 



where 



5T/T = 61+62 (2.12) 
0i = T-±[(P'Fj e -(P'F ti ) \e i , 

e 2 = \j S dXP"(r 1 )F ji e i . (2.13) 

Oi and 02 represent the intrinsic and Sachs- Wolfe effects, respectively. The latter can be divided into the 
ordinary Sachs- Wolfe effect Q S ac and the Integrated Sachs- Wolfe effect Qi sw > where 

®sac = \[{P"F) e -{P"F) ], 

6 Jst0 = ~ f * dXP'"{r,)F. (2.14) 
The second-order temperature perturbation is 

ST IT = J 1 (A e )[iji(A e )-r| - + r,e''] / " dXA^ 

2 lz J J\„ 

v{^ (2) ' + AWaW - A«" J* dXA^CX)} + ^d W \ (2.15) 



dX\ 



where (rj, x l ) — (A, A c — A) in the integrands and 

h(Xe) = ~\f ' dXP'F^e^, 
A (1) = ~PFi,-eV, 



2 

A (2) = -^[PLi + P 2 Mi + QN ij + C J i ]e i e i . (2.16) 



These expressions were derived in §3 of a previous paper. Q 



III. POWER SPECTRA OF SECOND-ORDER CMB ANISOTROPIES 

We consider primordial scalar perturbations with F defined by 

F = J dka(k)e ,;kx , (3.1) 

where the spatial average for a(k) is given by 

(a(k)a(k)) = (27r)- 2 ^(k),5(k + k), (3.2) 

with 

7V(k) - VFok-^k/k^-^ik), (3.3) 

where T s (k) is the matter transfer function [221 ] and Vfq is the normalization constant. 
Then the first-order temperature perturbations are 

ST/T = 6 P = -~ J dka(k) J ' dAF'(?7)(M 2 e lkx , (3.4) 

where x = re, kx = fcr/i, /i = cosdk and Ok is the angle between the wave vector k z and a unit vector e\ This 
equation can be rewritten as 

Q P = j dka(k){-i[(P'0 o + ^')o J Pi(/i)] + iEH)'( 2; + 1 ) e ^)^(^}, (3-5) 
J i 

where p (/i) is the Legendre polynomial and 

9 PW = ° d\P"'(r ] )j l (kr) - {k(P%(2l + + l)j l+1 (kr e ) - lji-i(kr e )]+ji(kr e ){P") e }. (3.6) 

In these equations, we have 7/ = A and r = A D — A. In the derivation of Eq. (|3.6[) , we used the relations [23l. l2~j] 

e ,:kx = e ifcr M = ^(_i)i(2Z + l)j,(Ar)fl(M) (3-7) 

and 

(2Z + l)/xP ; ( M ) = (/ + l)P m ( M ) + ZP i _ 1 ( M ). (3.8) 

The components of the unit vector e % are expressed as e 1 = sin 9 cos 0, e 2 = sin sin 0, e 3 = cos (9 with respect 
to x , x 2 , x axes, respectively. 

In order to derive the power spectra, we take the statistical average () for the primordial perturbations, and 
((5T/T) 2 ) is expressed for the first-order anisotropies as 

((ST/T) 2 ) = <(6p) 2 ) = (T )- 2 £^r^- ( 3 - 9 ) 

The power spectra C; are 

C*/ = (T ) 2 J dkk 2 V F {k)\U ( j}{k)\ 2 , (3.10) 
where To is the present CMB temperature and 

nf(k) = -(P") - k(P') eJ1 (kr e ) - (P") e j (kr e ) + J ° d\P"'j (kr), 

n ( p(k) = \k[{P') - {P') e ]3^{kr e ) - (P") e h(kr e ) + [ ' dXP'^^kr). (3.11) 
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For I > 2, we have 

H^(k) = fc(P')ei/ (1) (^e) - (P'OejH^e) + ^ dXP'" n(kr). 

In the derivation of C;, we used Eq. (|3.2p for (a(k)a(k)) and a mathematical formula for P(/i): 

Pi{[i)Pi<{[i)dn = 2/(21 + 1), for Z = l', I ± V , 

respectively. For the two directions with unit vectors ei and e2, we have the correlation 

27-1-1 

(r ) 2 (6 P (e 1 )ep(e 2 )) = V —^—dP^cos/3). 

z — ' 47T 

i 

where the product eie2 is equal to cos/3. 

For the second-order temperature anisotropies, we obtain from Eqs. (|2.15p . (|2.16[) and (|3.1[) 



(3.12) 



(3.13) 



(3.14) 



T/T = J J dkdka(k)a(£)|i J "J ' dXdXP' (X)P' (X)(k e k e ) 2 e t ^ x+ ^ 



dX 



ip'(A)(3fc e fc e + (k e ) 2 + (he) 2 - ^kk 



+ -Lp(A)P'(A) (l9/« e /<: e kk-14(fc e /<: e ) 2 -6(fc e A:) 2 -6(fc e /<:) 2 -3(kk) 2 



3£ 2 (fc) 2 



112 
1 
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— Q\x)(k e + k e f - (kk) 2 /(k + k) 



i(k+k)x 



r/A 



P"(A) / dXP(X)(k e k e f 



5 i(kx+kx) 



(3.15) 



where k e stands for ke, P'(A) = dP(X)/dX, and P'(A) = dP(X)/dX. Here r and A^' in Eq. ((2~l"5j) were 
neglected, because we pay attentions to the Sachs- Wolfe effect after the recombination epoch, and also the term 
with C\ in Eq. (|2.16p was neglected, because the contribution of gravitational radiation is very small. It is found 
from the above equation that the average (82 T/T) does not vanish in contrast to the vanishing first-order one 
({Si T/T}), and we obtain 



(ST/T) = I I dk(2ir)- 2 P F (k) 



dXdXP'(X)P'(X)(k e fe 



4 ik(x-x) 



dX 



\o J A 
1 



P'(A) (-(fc e ) 2 + - fc 2 ) + P(A)P'(A) (k 2 - 2(k e ) 2 ) (k e f 



+ 2 f ° dXP"(X) I dXP(X)(k e ) 4 e lk ^-^Y 



(3.16) 

where k e = ke = kfi and dk — dkk 2 d<pd^L. Here we have the relations x = re, x = fe, and so kx = k e r = 
kr(i, — kx = fcr(— /i). Thus e jk ( x ~ x ) is expanded as 

e ,k(x-x) = J2(-i) l (2l + l) Jl (kr)P l ( i i) x^H)''(2Z' + l)jV(fcf)fi'(-/i) 
l l' 
= £(-*)' **' {2l + l)(2l' + 1) ji(kr)ji>(kf) Pi(p)Pv(p). (3.17) 

i,l> 

Here for the reduction of Eq. (|3 . 16|) , we derive the following relations using Eq. (|3 . 8[) 

J> MP = £i, (1) P, 5> m 2 P = £ and ^ 3 P = £i, (3) P, (3.18) 



where 



^ (1) (fcr) 



Z 

2Z - 1 



ji-i(kr) 



Z + l 
2Z + 3 



3i+i(kr), 



G 



3?\kr) = 

+ 
+ 

Then we have 



(I - 1)1 



ji-2(kr) 



1 



(2/ — 3)(2Z — 1) 

(/ - 2)(l - 1)1 
(2/-5)(2/-3)(2/-l) 

I 2 (l + l) 2 



21 + 1 

jis(kr) H 



/ 2 



(« + l) 2 



2/- 1 (21 + 3) 



ji(kr) 



(Z + l)(/ + 2) . 



(2/ + 3)(2Z + 5)- 



I 



21 - 1 



{I -If 



(21 - 3) (21 - 1) 



(2Z-1)(2Z+1) (2Z + 1)(2Z + 3) 
(l + l) 2 (1 + 2) 2 



(2Z + 1)(2Z + 3) (2Z + 3)(2Z + 5) 



3i-i(kr) 
ji+i(kr) 



l + l 



2Z + 3L(2Z-1)(2Z+1) 

(I + l)(l + 2)(l + 3) 
(21 + 3)(2l + 5)(2l + 7) 



ji+3(kr). 



M 4 e ik (x -x) = J2(-i) l i 1 ' (21 + l)(2l' + l)jf ) (fcr)j i ( , 2) (fcf)P i (M)^'(M). 
Using these relations and executing the integrations in Eq. (|3.16p with respect to 4> and we obtain 
(ST/T) = (27r)- i y"dfcfc 2 ^(k){ifc 4 ^(2? + l)[^ Ae dAF\A)jf ) (fcr)] 2 

+ — (P e -Po)k 2 + — (P 2 ~P 2 )k 4 

24 v ; 840 e ° 



(3.19) 



(3.20) 



^ 4 ^(2Z + 1) / ' d\P"(X)j\ 2 \kr) f dXP(X)jl 2) (kf)}, 

1 " An J Xn 



(3.21) 



where we used the relation (|3.13[) . and P e — P(X e ) and P Q = P(X a ). 

When we take into account also the third-order anisotropies T/T, we have 

ST/T = ST/T + ST/T + ST/T 

= ST/T+ (S T/T) + Q PP + S T/T, 

12 3 

where O pp = Si T/T - (S 2 T/T). Then the total average of (ST/T) 2 is expressed as 

((ST/T) 2 ) = (Bl) + ({ST/T)) 2 + (B 2 ) + 2(5T/TST/T). 

F 2 FF 1 3 



(3.22) 



(3.23) 



The last term is of the same order as the second and third terms. At present we have not obtained any concrete 
expression of third-order metric perturbations yet and so the last term is not analyzed, while the formula for 
S3 T/T and formal solutions have recently been derived by D'Amico et al.[l6j and the perturbative equations to 
the third order were studied by Hwang and Noh[l7| with respect to the relativistic-Newtonian correspondence. 
Here (S2 T/T) is the monopole component without angular dependence and Q pp is the renormalized second-order 
temperature fluctuation, which have been discussed by Munshi et al.|| 
Now let us make a reduction of (Q pp )- It is expressed as 



(e 2 pp ) = ((jT/T(a(k)a(k)) - (5 2 T/T))(S 2 T/T(a(k)a(k)) (ST/T))), 



(3.24) 



where & T/T(a(k)a(k)) is given by Eq. (|3~T5| and <fe T/T(a(k)a(k)) is obtained from Eq. (|3~T5|) by re_placing 

k, k by k, k. The averaging process in Eq. J3.24|) is performed in the two sets: (a(k)a(k)), (a(k)a(k)) and 

(a(k)a(k)}, (a(k)a(k)}, and the average process for (a(k)a(k)) and (a(k)a(k)) is excluded by subtracting 

(S2T/T) from S2 T/T(a(k)a(k)) and S2 T/T(a(k)a(k)). When we sum these average values in the above two 
sets and symmetrize the expression with respect to k and k, we obtain from Eq. (|3.24l) 



dkdk(kk) 2 d(/) k d(l) k diidfi(2iT)- 4 P F (k)VF(k) 



{If J ' dXdXP'(X)P'(X)(k e k e ) 2 e l ^+^ + J ^ dx\^P\X)(3k e k e + (k e ) 2 + (k e f -\kk 
^P(A)P'(A)(l9fc e fc e kk - U(k e k e ) 2 - 6(k e k) 2 - 6(k e k) 2 - 3(kk) 2 + 3k 2 (kf 



i(k+k)> 
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-kk 

2 



~y " d\P"(X)J P(X){k e k e ) 2 e l ^+^} 

[\J ° J ° d V dfjP'(r 1 )P'(f])(k e ~k e ) 2 e^y+^ + J ° d V [^P'( V )(3k e k e + (k e ) 2 + (fc e ) 2 - 1 
^P(v)P'(v) (l9k e k e kk - U{k e k e ) 2 - 6(kk) 2 - 6(k e k) 2 - 3(kk) 2 + 3fc 2 (fc) 2 ) j e~ j(k+£)y 
~ / ^"(77) P dfjP(fj)(k e k e ) 2 e' l ^+^}, (3.25) 

" Ao " An 



where x, x and y, y are functions of A, A and 77,77, respectively, and we neglected the terms with Q and N'J, 

because we have Q/P 2 < 10~ 2 always and they are very small. Especially Q vanishes in the case A = 0. Here, 
when we consider an orthonormal triad vector e^, e}oy and (— e l ), the components of k and k with respect 

to this triad are expressed as k = k(sin9k cos fa, sinOk sinfa,cos9k) and k = k{svn9k cos fa, sin 6k sin fa , cos 9k ) , 
and so [i — cos 9k and /2 = cos 9k ■ 

Next let us take a notice of terms with kk and (kk) 2 which can be expressed as 

kk = kk [sin 9k sin 9k cos( fa — fa) + cos 9k cos 9k], 
(kk) 2 = (fcfc) 2 {(cos^cos^.) 2 + i(sin0 fc sin0 fe ) 2 [l + cos2(0 fe -0 fe )] 

+ 2 sin^fe sin^fc cos6*fc cos^fe cos(0fe — 0fc)}. (3.26) 

By integrations with respect to fa and fa,_we obtain J £ dfadfa cos(fa — fa) — J J dfadfa cos 2(0^ — fa) = 0, 
while / / fafa[cos{fa - fa)} 2 = f f d fadfa [cos 2(fa - fa)} 2 = (2tt) 2 /2. 

Executing integrations in Eq. (|3.25[) with respect to fa and (f>k , and using the relations Eq. (|3.7[) and 

e - iky = e -ik Sf * = J2(-i) m (2m + l)j m (ks)P m (-n) = ]Ti m (2m + l)j m (fcs)P m ( M ), (3.27) 

m tci 

we can, therefore, reexpress Eq. (|3.25[) as 

((©pp) 2 ) = E E E E i^Ml'mm' + Ml'mm' + 4w) ' (3.28) 

I V m m 1 

where the terms A I Wmm ,, A\l mm i and , come from the terms without fa and fa, the coefficients of 

cos 2 (0fe — fa), and the coefficients of cos 2 2(fa — fa), and their lengthy expressions are shown in Appendix A. 

Moreover let us replace the terms of [ijiPi, li 2 jiPi and fi 3 jiPi by jl Pi, j\ Pi and Pi using Eqs. (|3 . 18|) 
and (|3.19p and execute integrations with respect to fj, and p, using Eq. (|3.13|) . Then we obtain 

<( pp) 2 > = E E + B «' + Bu' 1 ) > ( 3 - 29 ) 



where 

= |(2tt)- 2 (2Z + 1)(2Z' + 1) /" J dkdk{kk) 2 V F (k)P F (k) 

^ j K dXd\P\\)P l {\)j\ 2) {kr)j[?\kf) + 2{kk) 2 f ' dXP"(X) [ X 

A„ JA iA JA„ 



x \{kkf / / dXdXP\X)P\X)j\ 2 \kr) 3 f{kf) + 2{kk) 2 dXP"(X) dXP(X)j (2) (kr)^ 2) (kr) 
Ja„ ja Ja Ja, 

^'(^(^fcjf^Mj'^C^) + fc'jf \kr)j v {kr) + k 2 ji(kr)jl 2) Ckr)) 



+ ^P(X)P'(X)(kk) 2 (ji(kr) Jl ,(kr) - Zj\ 2 \kr)j v (kr) - 3j,(fcr)j, ( , 2) (kr) 
+ 



ltf\kr)jP (kr))]} 2 , (3.30) 



= ^( 27r ^ 2 ( 2/ + 1 )( 2 ^' + 1 ) / I dkdk(kk) 4 V F (l;)P, (I;) 
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d\P' 



'(A)[(j,(*r) -tf\kr)) - ^ k kP(\)(jP(kr) - j?\kr))]j v (kr) 
d V P'(r,)[(j v (k S )-jP(ks)) - ^kkP( V )(jlPCks)~J^Cks))]jiCks), 



(3.31) 



P.]} 1 = 2 (^) 2 ( 27r ) _2 (^ + 1 )(^' + 1 ) / ( dkdk{kkfV F {k)V F {k) 



Ac 



dAP(A)P'(A)[ii(fcr)j I ,(fcr) - jf \kr)j v {kr) - ji{kr)j\?\kr) + j{ 2) (kr)tf ] {kr)\ (3.32) 



where r = A G — A, f = A D — A, s = rj — rj, r/ e = A e , and i] — A G . 

Next let us consider two directions with unit directional vectors ei and ei. If we define /^i,/x 2 and (3 as 
fii = kei/fe, H2 = ke 2 /fc and cos/3 = e!e 2 , respectively, we have a mathematical relation 



d(j) k dO k sin k Pi{m)Pi{n 2 ) 



47T 

27+1 



(cos (3). 



Using Eq. (|3.33p it is found that the correlation between 8 pp 's in the two directions is expressed as 
<e w ( ei )e pp (e 2 )} = \Y,H{ B tv +B $ + S" 7 )^(cos/3)PKcos/3). 



/ v 



Here let us expand PiPy by P n as 



Pt(cos/3)P^(cos/3) = J] Wn(cos/3). 



(3.33) 



(3.34) 



(3.35) 



The derivation of the coefficient bw n is shown in Appendix B. Then the correlation reduces to 

(r ) 2 (8 pp (e 1 )e pp (e 2 )) = £ ^±!c( 2 >P„(cos/?), 



(3.36) 



where 



C (2 



4tt 



2n 



T (T ) 2 ^(p/ ; , + P/ ( { + P^)^ n . 



(3.37) 



The expressions for ($ 2 T/T), ((9 PP ) 2 ) and (e pp (ei)6 pp (e 2 )} in Eqs.dSHH), ([3T2^]) and ([336]) are our new result 
which will be useful to derive the second-order power spectra. 



IV. CONCLUDING REMARKS 



(2) 

In this paper we derived the average value [S2 T/T) and the second-order power spectra C„ of CMB 
anisotropics due to primordial random density perturbations, which include two random variables a(k) and 
a(k), using the average values of the products of a(k). The average value (S2T/T) does not vanish and 

has no angular dependence. It should be regarded as the monopole component of temperature fluctuations 

(2) 

and the second-order angular correlation is described by the power spectra C« . Since we have not derived 
(<5i T/T S3T/T) yet, our analysis of second-order power spectra is incomplete. But we think (6 pp ) may repre- 
sent the essential feature of second-order power spectra, that is, their Z-dependence. In the next step we will 
analyze {Si T/T S3T/T) with the third-order metric perturbations for completeness. 

In the case A = 0, the first-order temperature fluctuations due to the integrated Sachs- Wolfe effect vanish, 
while the second-order ones do not vanish and play a dominant role in the i nteg rated Sachs- Wolfe effect, which 
should not be neglected. In the cosmological local void model f LVM) fisl ITg . I20I [2l| . the exterior region is 
described in terms of the Einstein-de Sitter model, and so the main part of the integrated Sachs- Wolfe effect in 
LVM is of second-order. In the interior region we use open low-density models with fio < 1, in which we have 
nonzero first-order Sachs- Wolfe effect. It is important for showing the observational reality of LVM to derive 
them. 
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In the case A ^ 0, the first-order temperature fluctuations due to the integrated Sachs- Wolfe effect decrease 
rapidly with the increase of the redshift z, but the second-order temperature fluctuations due to the integrated 
Sach- Wolfe effect decrease more slowly, and so the latter fluctuations may be dominant over the first-order 
fluctuations at the early stage. This situation is explained in a separate paper[25j using a simple model of 
density perturbations. Thus a characteristic behavior of CMB power spectra due to second-order temperature 
fluctuations will be measured through the future precise observation and bring useful informations on the 
structure and evolution of our universe in the future. 
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APPENDIX A: FIRST EXPRESSIONS OF A\ v mm , , A\l mm , AND A^ mm , 



A\ Vmm , = ^(2^)- 2 (-^) (^+^ ' ) ^ (m+ " l ' ) (2^ + l)(2r + l)(2m+l)(2m' + l)y J dkdk(kk) 2 d^dp,V F (k)V F (k) 
x [ "dx\ [ " dAP'(A)P'(A)(fcfe) 2 (///i) 2 ^ 



x ^J\ 

+ l -P{\)P'{\){kkf{2(^f - 6( M 2 + /2 2 ) + 3 - |(1 - M 2 )(l - m 2 ))~ 
x j l (kr)j v (kr)Pi(fi)P v (fi)+2P"(X) j d*J>(X)(kkiifi) a jt(kr)j v (fef)P,( M )Pj'(/2)} 



"2 



x p dr,{ J"" d77P'(? ? )P'(7y)(/cfc) 2 (^) 2 j m (fc S )j m ,(fcs)F m ( A/ )P m -(^) + J>fa)(|fcfc/i/i + fc 2 /x 2 + Pp, 2 ) 
+ ±P( V )P'( V )(kk) 2 (2( M/ 2) 2 - 6(/i 2 + /I 2 ) + 3 - |(1 - fi 2 )(l - m 2 ))" 

x j m (ks)j m *(ks)P m (Li)P m -(p,)+2P' , (ri) df}P(f})(kkw) 2 j m {ks)j m ,(ks)P m ([i)P m ,{fi)y (Al) 

Jrjo 

Ml> mm > = ^{2Tr)- 2 (-i) (l+l '\( m+m '\2l + l)(2l' + l)(2m+l)(2m' + 1) J J dkdk(kk) 4 dfidp,V F (k)V F {k) 

f x ° r 26 i 

x / dAP'(A)[l-yP(A)fcfc^J^(fcr)>(fcr)(l- A1 2 )P i ( A1 )P r (/2) 

f' 1 " r 26 - 1 

x / drjP\ri) 1-—P(r])kk^ j m {ks)j mr (ks)(l-ft 2 )P m (fx)P m '(fl), (A2) 
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Afl'Lm' = 2 (yu) ( 2 ^) _2 (- J ) (W ' ) * (m+m ' ) ( 2Z + 1 )( 2r + 1 )( 2m + 1 )( 2TO ' + 1 ) j j dkdk(kk fdfldftP F (k)P F (k) 
x [' dXP{X)P'{X)j l {kr)j v {kr){l - M 2 )(l - m 2 )P ; (^)Pp(m) 



A 



X 



>lo 



driP(r))P'{r))j m (ks)j m ,(ks)(l - M 2 )(l - Jl 2 )P m {n)P m '{p), (A3) 



where r = |x|, s = |y|, r = X a - X, r = X — X, s = rj — T), s = r/ Q — r), r\ e = X e , and r\ = X a . 



APPENDIX B: DERIVATION OF b„, 



11' n 



Since the Legendre function P m (z) is a finite power series of z 3 with integers j(< m), P m (z)P m '(z) can be 
also expressed as a finite power series of z J with integers r satisfying < j < m + m'. On the other hand, z n 
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is expressed as a finite series of P m {z) with < m < n, so that P m (z)P m '(z) can be expressed as a finite series 
of P n (z) with < n < m + m! . Their examples are 

(P 2 ) 2 = JPo + ^ 2 + gp. 

P1P3 = yP 2 + yP 4 , 

PP - 9 P + 4 P + 10 P 
P 2 P 3 - 35 P i+i5 P 3+21 P 5. 

In these cases the coefficients in Eq. (|3.35j) are 



6121 


= 2/5, 


6123 = 


3/5, 




&220 


= 1/5, 


6222 = 


2/7, 


6224 = 18/35, 


&132 


= 3/7, 


^134 = 


4/7, 




&231 


= 9/35, 


6233 = 


= 4/15, 


6235 = 10/21 



The general formula of these series are classified into the following three cases. 
(1) even- even type 



(Bl) 



(B2) 



m+m 



PlmPlm' — E] b 



1m 1m' Is ^1s, 



s=0 



Jim 1m' Is 



4s+ 1 
o2m+2m 



7 E £(-l) m+m ' +J ' +j '(2j + 2/)!(2m + 2j)!(2m' + 2j')\ 



3=0 j=o 

I (m + i)!(m - j)!(2j)!(m' + /)!(m' - /)!(2j')!(2j + 2j' - 2s)!!(2j + 2j' + 2s + l)!!l , (B3) 



(2) odd-odd type 



m+m +1 

Plm+\Plm' + 1 = &2ro+l 2m' + l 2s P 2 S) 



>2m+l 2m' + l 2s 



4s + 1 

22m+2m'+2 



E J2(-l) m+m ' +j+j '(2j + 2f + 2)\(2m + 2j + 2)!(2m' + 2/ + 2)! 



(3) even-odd type 



J'=0j'=0 

+ l)!(m - i)! (2j + l)!(m' + f + \)\{m' - j')\{2j' + 1)! 
x (2j + 2f -2s + 2)!!(2j + 2/ + 2s + 3)!!] , 

m+m 



(B4) 



P2mP2m' + l — E; 1m' + 1 2s+l P 2 s+1, 



s=0 

4s + 3 



32m Im' + l 2s+l 



22m+2m 



— EE (-l) OT+m ' +i+J "(2j + 2j' + l)!(2m + 2j)!(2m' + 2f + 2)! 



3=0 j'=0 

/ (m + i)!(m - j)!(2j)!(m' + j' + l)!(m' - j')\(2f + 1)! 
x (2j + 2/-2a)!!(2j+2/ + 2a + 3)! 



(B5) 



where the rules of ! and !! are 



n! = 1 ■ 2 • • • • (n — 1) • n for n > 1, 
0! = 1, 

nl — oo for n < 0, 



(B6) 



and 



(2n)!! = 2-4----(2n-2)-2n for n > 1, 
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(2n 



1)!! = 1-3- 
0!! = 1, 
nil = oo 



(2n 



for 



- 1) -2n+ 1 
n < 0. 



for 



n > 0, 



(B7) 
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